CENTER VORTICES, MAGNETIC CONDENSATE 
AND CONFINEMENT IN A SIMPLE GAUGE SYSTEM 



F. GLIOZZI AND M. PANERO 

Dipartimento di Fisica Teorica deU'Universitd di Torino, 
Istituto Nazionale di Fisica Nucleare, Sezione di Torino, 
via P. Giuria 1, 1 10125 Torino, Italy 

AND 

P. PROVERO 

Dipartimento di Scienze e Tecnologie Avanzate 
Universitd del Piemonte Orientale, Alessandria, 
Istituto Nazionale di Fisica Nucleare, 

gruppo coHegato di Alessandria, 1 15100 Alessandria, Italy 



Abstract. The confining mechanisms of 't Hooft and Mandelstam have a 
simple microscopic realization in 3D Z2 gauge theory: the center vortex and 
the magnetic monopole condensation are associated, in the set of config- 
urations contributing to the confining phase, to the presence of two kinds 
of infinite clusters. These generate the area law of the large Wilson loops 
and the universal finite size effects produced by the quantum fluctuations 
of the bosonic string describing the infrared behavior of the flux tube. 



1. Introduction 

Two kinds of topological excitations were proposed by 't Hooft as possible 
composite fields condensing in the disordered phase of gauge systems to 
produce confinement: magnetic monopoles [1] and center vortices [2]. 

Condensation of magnetic monopoles implies dual superconductivity of 
the vacuum and gives the well-known physical picture of confinement in 
terms of dual Abrikosov vortices which describe the flux tubes. 

Center vortices are string-like structures which are created by gauge 
transformations with a non-trivial homotopy. Their condensation produces 
a very efficient disordering of the gauge configurations. 
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In the last few years, many lattice studies have been made on these two 
subjects in the last few years (see for instance the two review talks of Di 
Giacomo [3] and Olejnik [4] at this Conference). 

These two kinds of confining mechanisms have a particularly simple 
implementation in the 3D Z2 gauge model, where it is possible to study 
this phenomenon not only in terms of the effective degrees of freedom, 
using for instance the underlying string description of the infrared behavior 
of the flux tube, but also at a microscopic level, i.e. in terms of Z2 lattice 
configurations: it is possible to pick out the degrees of freedom which are 
relevant for the infrared dynamics of the confining phase. 

It turns out that the condensation of the above-mentioned topological 
excitations manifests itself at the microscopic level with the appearance 
of two (complementary) infinite clusters of links in the dual lattice. Both 
produce not only area law decay of the vacTium expectation value of large 
Wilson loops, but also the universal finite size effects which are predicted 
by the free bosonic string picture of the flux tube. 

2. Monopole condensation 

The three-dimensional Z2 gauge model on a cubic lattice is deflned by the 
action 

Sif3) = -f3Y,Un ,Ua = llUe (1) 

o tea 

where [/^ = ±1 is the Z2 gauge fleld associated to the links i of the lattice 
A. The magnetic monopoles live in the dual lattice A. In order to create 
a monopole in a site x G A corresponding to the center of an elementary 
cube of A it is sufficient to draw an arbitrary, continuous line j{x) joining 
a: to 00 and flipping the sign of the coupling of the plaquettes crossed by 
7. This flipping is generated by the non-local operator 

*^(x)=exp(-2/3^i7a) . (2) 

□£7 

As a consequence, the flux across any closed surface around x is equal to -1: 
this monopole fleld ^'y{x) creates one unity of Z2 flux. Thus the monopole 
is associated to a point-like topological defect in the ensemble of the gauge 
conflgurations. 

Monopole condensation occurs when ^' acquires a vacuum expectation 
value different from zero: 

monopole condensation <^=^ {^) 7^ . 

Why this condensation should imply conflnement? A useful piece of in- 
formation comes from the Kramers- Wannier duality. This transformation 
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maps the Z2 gauge theory into the Ising model described by the Hamilto- 
nian 

H{P) = -^J^aiay ,~^ = -\ log tanh (3 (3) 

{iy) 

with {xy) ranging over the hnks of A and the spin variable Ux taking the 
values ±1. As a consequence one can easily show that 

= {ax)lsing (4) 

Thus the magnetic condensation is associated with the spontaneous Z2 
symmetry breaking of the Ising model. 

The same transformation maps the Wilson loop W{C) associated to 
any closed curve C G A into the corresponding 't Hooft loop W{C) of the 
Ising model: 

{W{C)) = {W{C))i,ing (5) 

with 

W{C) = exp I -2/3 OiUj | , OS = C (6) 

whCTe S is an arbitrary surface bounded by C. 

W{C) creates an elementary Z2 flux along C which manifests itself as 
a topological defect or a frustration in the ensemble of the Ising configu- 
rations: the product of the link variables along any loop C G A having an 
odd linking number with C is equal -1. 

Confinement requires that {W{C))ising should decay with an area law 
in the ordered phase of the Ising model. To understand this property at a 
microscopic level it is convenient to resort to the Fortuin Kasteleyn (FK) 
random cluster representation of the model: 

Z = ^e-^(^) = J]t;^«2^<' , (7) 

{a} GCA 

where v = e^^ — 1, the summation is over all spanning subgraphs GCA. 
bo is the number of links of G, called active bonds, and cq is the number 
of connected components, called FK clusters. The ordered phase is charac- 
terised by the presence of an infinite, percolating FK cluster. 

When there are frustrations in the system (for instance those generated 
by W{C)) the summation in Eq.7 is constrained. Only those FK clusters are 
allowed for which no closed path within the cluster is linked to an elemen- 
tary Z2 flux [9] . This is a microscopic realization of a sort of dual Meissner 
effect: the FK clusters behave like pieces of dual superconducting matter of 
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Figure 1. Density of the largest center vortex (squares) and FK (circles) clusters. The 
densities of the second largest clusters of the two types are also shown. 



type I and no Z2 flux can go through them. Introducing a projector wc{G) 
on the configurations G C A which takes the value whenever there is a 
FK cluster (at least) linked to C, and the value 1 in all the other cases, 
yields the very useful identity 

, number of allowed config. 

C Ismg total number of config. 

Assume that there are only FK clusters of finite size ( this is the case of the 
symmetric phase). If C is much larger than the mean size of the clusters, 
the configurations contributing to W{C), i.e. those with wc{G) = 1, are 
characterised by the fact that there is no cluster linked to C. The weight of 
this class of configurations, when compared with the total ensemble G C A 
is clearly suppressed by a factor e""^'''^'*, where p{C) is the length of C. 
Thus the Wilson loop obeys a perimeter law in this phase. Conversely we 
can argue that a decay of large Wilson loops with an area law implies the 
presence of an infinite FK cluster. Quite remarkably, this kind of reason- 
ing, based only on concepts well defined at the microscopic level, leads 
to the very important, unavoidable, conclusion that confinement requires 
monopole condensation. This is just the reverse of the usual 't Hooft con- 
finement criterion, which can be justified using plausibility arguments based 
on the effective description in terms of (dual) Abrikosov vortices. 

The above arguments lead us to conjecture that the infinite FK cluster 
encodes the most relevant part of the infrared dynamics of the confined 
phase. More precisely, starting from the the asymptotic form of the vacuum 
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expectation value of a large rectangular Wilson loop 



{W{R, T)) = exp {~aRT + p{R + T) + c) 



Vr 
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(9) 



r?(r) 



where ?7(t) is the Dedekind function and the term within square brackets 
is the universal contribution due to the quantum fluctuations of the flux 
tube, it is easy to prove a sort of no-renormalization theorem: the vacuum 
expectation value of W{R, T) in the ensemble of configurations obtained by 
those in thermal equilibrium erasing all the finite FK clusters has exactly 
the same functional form of Eq.9 and with the same value of a. The argu- 
ment goes as follows: let N be the total number of configurations and Ni 
those compatible with C. We have {W{C)) = j^. When the finite FK clus- 
ters are eliminated, the number N[ of compatible configurations increases 
(iV( > A'^i) because some graph with wciG) = is promoted to a compat- 
ible one. The perimeter law decay produced by the finite FK clusters can 
now be used to get A^i = e^"^'^'^-' A''(, which completes the proof. 

The transformation Ni — N[ reduces drastically the noise in the eval- 
uating of large Wilson loops and it has been used also in related models 
to study the universal shape effects related to the quantum fluctuations of 
the flux tube [5] . 

3. Clusters of center vortices 

Center vortices in Z2 gauge model arc constructed by assigning a vortex 
line in the dual lattice to each frustrated plaquette (i.e. Un = —1) in the 
direct lattice. Since the product of the six plaquettes forming a cube is 
constrained to be equal to 1, the resulting graph of center vortices (CV) in 
A has even coordination number. This has to be contrasted with the FK 
clusters of the previous section which are arbitrary subgraphs of A. 

The value of a Wilson loop W{C) in a given configuration is ±1 ac- 
cording to the number, modulo 2, of frustrated plaquettes of an arbitrary 
surface S bounded by the loop (9S = C), or, in the center vortex language, 
to the number, modulo 2, of vortex lines that are linked to the loop. The 
role of center vortices in confinement in this model is in a way trivial, since 
removal of all center vortices from Z2 gauge theory configurations simply 
removes all the dynamics by transforming every configuration into the triv- 
ial vacuum. Therefore center vortices are, in this sense, trivially responsible 
not only for confinement, but for the whole dynamics of the model. 

We are interested to a subtler issue, namely the effect of the finite 
CV clusters on confinement [6]. Also in this case it is easy to prove that 
confinement requires an infinite CV vortex, however it is no longer possible 



6 



to prove the no-renormalization theorem of the FK formalism, hence there 
is no reason to beheve that the string tension does not change under the 
transformation erasing the finite CV clusters. 

It is straightforward to verify numerically that in the confining phase of 
Z2 gauge system there is no ambiguity in finding a cluster of center vortices 
whose size scales linearly with the lattice volume for large enough lattices 
(sec Fig.l), while in the deconfined phase we found that the density of the 
largest cluster decreases rapidly with the volume. 

It has to be noted that the presence of the infinite cluster does not nec- 
essarily imply a percolation property of the central vortices. Confinement 
requires merely the presence of an infinite cluster. In the FK representation 
this requirement is sufficient to assure percolation of the infinite FK cluster, 
while there are regions in the confining phase where the infinite CV cluster 
does not percolate, perhaps because of the even coordination constraint. 

To study more carefully the relationship between the presence of the 
infinite cluster and the value string tension, we chose [10] to simulate the 
model at /3 = 0.74883, which is well inside the scaling region, and for which 
the value of the string tension is known with high precision from simulations 
of the dual model, that is the 3D (spin) Ising model [7,8] a = 0.01473(10) 

First, we verified the existence of an "infinite" cluster of vortex lines: 
for each configuration, we selected the largest connected component of the 
graph defined by the center vortices, and verified that the size of such 
component grows linearly with the volume of the lattice. 

The results of this analysis are shown in Fig. 1, where the size of the 
largest CV cluster divided by the lattice volume is shown to approach a 
constant for large lattices. The size of the second largest cluster is also 
shown: its size relative to the lattice volume tends to zero and the identi- 
fication of the "infinite" cluster is unambiguous. For comparison, also the 
corresponding FK cluster densities are reported. 

To test the relevance of the largest cluster of center vortices, we pro- 
ceeded as follows: first, we modified each configuration in the Monte Carlo 
ensemble by eliminating all the vortices not belonging to the largest cluster, 
and second, by eliminating, instead, the largest cluster only. The qualita- 
tive picture described above suggests that a non-zero string tension will be 
found in the first case but not in the second. 

An efficient method to extract the string tension from Wilson loop data 
generated by Monte Carlo simulations, which takes into account the string 
fiuctuation contribution, was introduced in Ref. [11]. One defines the ratio of 
the expectation values of rectangular Wilson loops with the same perimeter: 



(10) 



7 



Full configurations 




Figure 2. (a) ac:ft{R,n) as a function of 7? for n = 2,3,4,5,6 for the original config- 
urations, (b) Th same quantity for the configurations where only the largest cluster of 
center vortices has been left, (c) The effect of removing only the largest cluster. 



Using the asymptotic expansion defined in Eq.9 we get 

r{R,n) exp{an'^) F{n/R) , (11) 



where F{t) = ^Jr]{i)^/l — t/rj ^ij^^, so that one can define the following 
quantity 

which approaches the string tension for large enough R. The results of our 
simulations are reported in Fig. 2 

In conclusion, our results confirm the picture of confinement as due to 
the existence of an infinite cluster of center vortices: our choice of the Z2 
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gauge theory allows us to bypass all the problems related to the gauge- 
fixing and center projection that one encounters when studying the same 
issue in SU (N) gauge theories. Two important new facts emerge from our 
study: 

While the largest center vortex is responsible for confinement, since its 
removal from the configurations makes the string tension vanish, the string 
tension measured from configurations in which all the other clusters have 
been removed does not reproduce the full string tension of the original 
theory. Therefore small clusters of vortices, while unable by themselves to 
disorder the system enough to produce confinement, do give a finite con- 
tribution to the string tension of the full theory. This has to be contrasted 
with the FK cluster formulation, where the infinite cluster account for the 
whole string tension, because of a no-renormalization theorem. 

The quantum fluctuations of the flux tube survive the elimination of 
the small CV or FK clusters: the Wilson loop after deletion of all the small 
clusters show the same shape dependence as the ones of the full theory, 
which can be explained as originating by the fluctuations of a free bosonic 
string. 
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